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1. Introduction 



Confinement of color is a well established property of strongly interacting matter at low tem- 
peratures and densities, even if not yet fully understood starting from the QCD first principles. The 
existence of a high temperature phase transition to a deconfined state of matter has been predicted 
by numerical simulations of lattice QCD and is presently under investigation in Heavy Ion exper- 
iments. In the present study we address the question regarding the fate of confining properties as 
the low temperature - high density phase transition is crossed: that is relevant in order to charac- 
terize the nature of matter in compact astrophysical objects and more in general to understand how 
deconfinement at high densities compares to what happens at high temperatures. 

Numerical studies of QCD at finite density are notoriously difficult because of the sign prob- 
lem, which makes usual importance sampling simulations unfeasible: for that reason we have stud- 
ied the theory with 2 colors, where that problem is absent. No sensible differences are expected for 
the confining properties of the theory when going from N c = 2 to N c = 3, where N c is the number 
of colors: for that reason we believe that our study could be relevant also for real QCD. 

Indications about deconfinement at high density obtained so far have been based on the anal- 
ysis of the Polyakov loop [j|], which however is not a true order parameter for confinement in 
presence of dynamical fermions. Different order parameters can be constructed in the framework 
of specific mechanisms for color confinement. One successful mechanism is that based on dual 
superconductivity of the QCD vacuum ^ ||, ||], which relates confinement to the spontaneous 
breaking of a dual magnetic symmetry induced by the condensation of magnetic monopoles: in 
that context a disorder parameter can be developed [Q, ^] which is the expectation value of a mag- 
netically charge operator, (^#) 1 , and which has been successfully tested as an order parameter for 
color confinement both in the quenched theory [J7|, |8|] and in presence of dynamical fermions |Q [1C| ]. 
Similar parameters have been constructed elsewhere [11, 12, |13|, 14]. 

Our plan is to study the behaviour of in the whole T — fx plane, in order to char acterize the 
confining properties of the various phases in the QCD phase diagram. In Section ^| we will review 
the definition of the disorder parameter and present our strategy for its numerical study; some 
preliminary results concerning the theory with 8 flavors of staggered quarks will be presented in 
Section ||. 



2. The disorder parameter {Jit) 

^(x,t) is defined in the continuum as the operator which creates a magnetic charge in x,t by 
shifting the quantum field by the monopole vector potential (y — x) 




(2.1) 



Its expectation value, when discretized on the lattice, appears in the following form 

{Jt)=Z/Z\ (2.2) 

'We change the usual notation for the disorder operator, (ju), in order to avoid confusion with the notation for the 
chemical potential. 
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Z is the usual QCD partition function 

Z = J (&U)detM(n)e- ps (2.3) 

where M is the fermionic matrix and S is the pure gauge action, while the partition function Z is 
obtained from Z by a change in the pure gauge action S —> S which adds the monopole field to the 
temporal plaquettes at timeslice t. 

Being expressed as the ratio of two different partition functions, the numerical study of 
is a highly non trivial task: while numerical methods have been recently developed which permit 



a direct determination of (^m) [15], we will not use them in the present study since they involve 
the combination of several different Monte Carlo simulations, which in presence of dynamical 
fermions could be unpractical. We will instead study, as usual, susceptibilities of the disorder 
parameter, from which the behaviour of (^#) at the phase transition can be inferred. 

For instance, being interested in (JK) as a function of j3, as for the ji. = phase transition, one 
usually measures [|], ^, [7|] 

P = ^^) = ^Z-^lnZ=(S) s - { S) § (2.4) 

where the subscript means the pure gauge action used for Monte Carlo sampling. The disorder 
parameter can be reconstructed from the susceptibility p , exploiting the fact that (M) = 1 at j8 = 

<^)(j8)=expQ%(j8')dj3') . (2.5) 

In particular p ~ in the confined phase means (^) ^ 0, a sharp negative peak of p at the phase 
transition implies a sudden drop of (^#) and p diverging in the thermodynamical limit in the 
deconnned phase means that (^#) is exactly zero beyond the phase transition. 

At finite temperature and density we are interested in studying the behaviour of (^#) in the 
two parameter space (j8,ju), where pL is the chemical potential. For that reason we introduce the 
new susceptibility 

<*> s £ XAJ{) = ^-^' oZ = {N ' )s - {N ' h <2 - 6) 

where Nf is the quark number operator. The dependence of (^#) on the chemical potential /I can 
then be reconstructed as follows: 



, i u) = (^)(^,0)exp^p D ( At / )d J u / 



so that, if the starting point at = is in the confined phase ((^#)(j8,0) ^ 0), the behaviour 
expected for Pz)(A t ) in correspondence of a possible finite density deconfinement transition will be 
the same showed by p across the finite temperature transition. 

3. Numerical results 

We have used staggered fermions corresponding to Nf = 8 degenerate continuum flavors, with 
bare mass am — 0.07. Different lattices have been considered, with a fixed temporal extent Lf — 6 
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Figure 1: Results obtained for po 



and a variable spatial size (only L s = 8, 16 so far) in order to make a finite size scaling analysis of 
the phase transition. The critical value of j3 at ju = is j3 c — 1.59 JT6]]: we have made simulations 
at different values of the chemical potential /I and at a fixed value of j8 = 1.5 < j8 c , i.e. starting from 
the confined phase at pL = 0. An exact HMC algorithm has been used and standard actions both 
in the gluonic and in the fermionic sector: we have used trajectories of unit length and a step size 
variable from 10~ 2 to 4 • 10~ 3 : the reduced time step was necessary in order to keep a reasonable 
acceptance rate around and above the phase transition, as a consequence of the presence of small 
eigenvalues of the fermion matrix. The same problem was at the origin of a severe slowing down 
around the critical value of pL, which made the availability of the recently installed apeNEXT facil- 
ity in Rome essential in order to carry out simulations on the larger lattice (L s = 16). Simulations 
on the smaller lattice L s = 8 have been carried out both on the APEmille facility in Pisa and on the 
PC farm of INFN in Genova, using there a code adapted from the publicly available version by the 
MILC collaboration. 

Our results for as a function of /I are shown in Fig. [j]. pr> shows a clear peak at a critical 
value of the chemical potential jx c ~ 0.3. The peak deepens as the lattice volume is increased, sug- 
gesting the presence of a true phase transition at which (./#) drops to zero and confinement (dual 
superconductivity) disappears. The position of the peak is coincident with that of other susceptibil- 
ities and we show in Fig. ||the case of the Polyakov loop: our data are still noisy and do not show 
any clear size dependence, which however is not expected since the Polyakov loop is not an order 
parameter for the phase transition. 

A more detailed analysis of the disorder parameter around the phase transition can be carried 
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Figure 2: Results obtained for the suscepibility of the Polyakov loop. 



out as follows. At fixed T we can assume the following finite size scaling behaviour for (^#) as a 
function of /x. 

from which it is easy to derive 

Pd =4 / X(m-m c )4 /v ) (3-D 

Our data show a nice scaling with v ~ 0.66, as shown in Fig. [|[ which is compatible with a 
second order phase transition in the universality class of the 3d Ising model. 

4. Conclusions 

We have made use of the disorder parameter detecting dual superconductivity of the 
vacuum to inspect the confining properties of QCD with 2 colors at finite temperature and density. 
The study has been carried our with Nf = 8 flavors of staggered fermions. In order to analyse 
as a function of the chemical potential \i we have introduced a new susceptibility = 
d/dji ln(^), showing the presence of a phase transition at finite density where confinement (dual 
superconductivity) disappears; the transition is coincident with those signalled by other quantities 
such as the Polyakov loop. A preliminary finite size scaling analysis of the disorder parameter is 
compatible with a phase transition in the universality class of the 3d Ising model. 

After these preliminary results, we plan in the future to make a more extensive study of the 
disorder parameter in order to characterize the whole phase diagram of the theory with two colors, 
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Figure 3: Finite size scaling analysis for po 



with a particular interest in the region of low temperatures and high densities. To that aim, a com- 
bined study of both susceptibilities, p and po, could be particularly useful: indeed the knowledge 
of Vln(^#), with V = (d/dfi, djd\l), in the whole j3 — \l plane, could give information not only 
on the location but also on the direction of the critical line, thus permitting a more careful study of 
the phase diagram. 
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